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The high-energy QCD factorization for Deep Inelastic Scattering and for proton-nucleus 
collisions using Wilson line formalism and factorization in rapidity is discussed. We show 
that in DIS the factorization in rapidity reduces to the fc-p-factorization when the 2-gluon 
approximation is applied, provided that the composite Wilson line operator is used in 
the high-energy Operator Product Expansion. We then show that the inclusive forward 
cross-section in proton-nucleus collisions factorizes in parton distribution functions, frag- 
mentation functions and dipole gluon distribution function at one-loop level. 
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1. Introduction 

In Quantum Chromodynamics, factorization is a fundamental concept. It allows the 
applicability of perturbative methods to the calculation of scattering amplitudes 
and cross sections of hadronic processes. Well known factorization schemes are the 
collinear factorization and the fcx-factorization. These two factorization schemes are 
relevant when the hadronic matter is probed at two different limit: the Bjorken limit 
and the Regge (high-energy or small-x) limit. The collinear factorization is applied 
in the Bjorken limit where the dynamics of the process is governed by incoherent 
interactions. On the other hand, the fox-factorization is relevant at high-energy 
(Regge limit) where the dynamics of the process is dominated, instead, by coherent 
interactions. 

At high-energy (Regge limit) the scattering amplitude is suitably factorized in 
rapidity-space (see Fig. ([!])). and the coefficient- functions and matrix elements of 
non-local operators of the Operator Product Expansion at high energy, contain per- 
turbative and non perturbative contributions. The non local operators are Wilson 
lines: infinite gauge link ordered along the straight line collinear to the particle's 
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Fig. 1. Expansion of the T-product of two electromagnetic currents in terms of Wilson-line op- 
erators. The blue dotted lines represent the Wilson line operators. 



Here rj is the particle's rapidity. The evolution in rapidity of these operators is 
known as the Balitsky-equatiorlH a non-linear evolution equation which generates a 
hierarchy of coupled equations known as the Balitsky- hierarchy^ (for a review seeO . 
In the Color G lass Condensate formalism it coincides with the JIMWLK evolution 
equatiorPHE In the large N c limit, instead, the Balitsky-equation decouples and 
is written in a clo sed form that is known, in DIS the Balitsky-Kovchegov 

(BK) equatiorW The linear version of the BK equation is the BFKL equation-^. 

2. NLO photon impact factor for DIS 

Let us now illustrate the logic of the OPE at high energies when applied to the 
T-product of two electromagnetic currents which will be relevant for DIS process. 
The technique we are using is the background field technique: the T-product of 
the two electromagnetic currents is considered in the background of gluon field. In 
the spectator frame the background field reduces to a shock wave (for review see 
[2]). In DIS, in the dipole model, the virtual photon which mediate the interactions 
between the lepton and the nucleon (or nucleus), splits into a quark anti-quark 
pair long before the interaction with the target. The propagation of the quark anti- 
quark pair in the background of a shock wave, reduces to two Wilson lines. If the 
quark fluctuate perturbatively in a quark and a gluon before interacting with the 
target, then the number of Wilson lines increases. Formally, we can write down the 
expansion of the T-product of two electromagnetic currents in the following way 



velocity 




(1) 
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where U x = Pexp(ig J dx + A _ (x + + xj_) is the Wilson line. In Eq. ([3]), the coefficient 
lltv represents the leading order impact factor, while the NLO impact factor is given 
by the coefficient i™ . In QCD, Feynman diagrams at tree level are conformal 
invariant. The LO impact factor is indeed conformal invariant and it can be written 
in terms of conformal vectors k = 55^(^7 — x 2 P2 + %±) — 3^ (*F — V 2 P2 +V±) an d 

^^;[2( K .Ci)(«-C 2 )-- 2 (Ci-C 2 )] + 0(a s ) (3) 

Although the NLO impact factor is also made by tree level diagrams, it is not 
conformal invariant due to the rapidity divergence present at this order. Since we 
regularize such divergence by rigid cut-off, we introduce terms which violate confor- 
mal invariance. In order to restore the symmetry we introduce counterterms which 
form the composite operator. The procedure of restoring the loss of conformal sym- 
metry due to the regularization of the rapidity divergence by rigid cut-off, is analog 
to the procedure of restoring gauge invariance by adding counterterms to local op- 
erator when the rigid cut-off is used instead of dimensional regularization, that 
automatically preserve gauge symmetry, to regulate ultraviolet divergence at one 
loop order. 

In Eq. ([3]), the composite operator is 

[Tr{^^}] C ° nf = Tr{UlU%} (4) 

V7Y J z 13 z 23 jv c ^13^23 

The parameter a is the analog of fip in the usual OPE. Note also that at this order 
the operator proportional to the NLO impact factor does not need to be modified. 
It would get a counterterm at NNLO accuracy. Using, then, the composite operator, 
the NLO impact factor is conformal invariant and it can be written entirely in terms 
of the conformal vectors we defined above. See Ref.[5]for its explicit expression. Such 
result is an analytic expression of the photon impact factor in coordinate space which 
is relevant for DIS off a large nucleus where the n on lin ear operator appearing at 
NLO level is relevant at high parton density regime^^. 



2.1. NLO photon impact factor for DIS in the fcx factorization 
scheme 

What we are interested in is the NLO impact factor for BFKL pomeron in mo- 
mentum space. Thus, our next step, before proceeding to the calculation of the 
Mellin representation which is a useful intermediate step to get the momentum 
representation, is to obtain the linearization of result in coordinate space in the 
non-linear case. Diagrammatically, the linearization of the non linear terms of the 
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Fig. 2. Diagrammatic representation of the 2-gluon approximation applied to Wilson line opera- 
tors proportional to LO and NLO impact factor respectively. 



non- linear equation is given in Fig. (|2.1[) . To get the right-hand-side (RHS) of Fig. 
(|2.ip we applied the 2-gluon approximation, thus, reproducing the diagrams of the 
NLO Impact factor of the usual perturbative QCD. It turns out that the coordinate 
representation of the NLO impact factor in the linea rized case can be written as a 
linear combination of five conformal tensor structure d^ H 

The projection of the impact factor on the Lipatov eigenfunctions with con- 
formal spin is related to the unpolarized structure function for DIS. While the 
projection on the Lipatov eigenfunction with conformal spin 2 is related to the po- 
larized structure function. The result of the Mellin representation can be found in 
Ref. 1121 Once we have performed the Mellin representation, we are ready to perform 
the Fourier transform in momentum space. The result of the Fourier transform iJ^l 



Vy ' ±; 32y 7ri/(l + i/2)cosh 2 7riAQ 2 / 

{[(j^X' + f + ^'K 



where 



pjiv 



1ix.1v , 



r 2 



q* v q-P2s \ q-P2 



(g^ ig* - p£-«-) (<f 1 - z^ 2 - ft-*-) 

q-P2 q-?2 
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U + *ir -P2 ){g +w -P2 J 

q-P2 q-P2 
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and 



= ^i(7) + 2^(2-7)-2^(4-27)-V(2 + 7), 




In order to obtain the full expression in momentum space of the NLO DIS 
amplitude, we need to perform also the Fourier transform of the NLO linearized BK 
equation for the dipole form of the unitegrated gluon distribution V(z) = z~ 2 U{z) 
where U(x,y) = 1 — NT 1 ti{U(x±U^(y±)}. The fcx factorized formula of the NLO 



where I PlU (q,k±) is given in Eq. ([5]) and the evolution of the operator V 0m= i fl (fci) 
can be found in Ref. [T2land it is in agreement with result of Ref . [TBI HH (see also Ref. 
H5l H6l HTf . The &t factorized formula (J6j) is a consequence of the factorization in ra- 
pidity of the DIS amplitude obtained through the high-energy OPE with composite 
Wilson line operator in the 2-gluon approximation. In other words, the factorization 
in rapidity of the DIS amplitude naturally reduces to the /cx-factorization when the 
2-gluon approximation is applied. 

3. One-loop factorization for inclusive hadron production in pA 



In proton-nucleus collision in the forward region the two factorization schemes, the 
collinear factorization and the rapidity factorization, enter on the same footing: the 
proton is treated as a diluted system which, using collinear factorization, emits a 
quark or a gluon that eventually scatters off a dense target like a large nucleus. At 
this point it is important to know whether the naive LO factorization formula of the 
cross-section written as a convolution of the parton distributions and fragmentation 
functions and dipole gluon distribution holds also at NLO level. 

To illustrate the one loop calculation, let us consider the quark-channel. The 
result of the calculation of the real and virtual Feynman diagrams at one loop order 
(shown in Fig. (J3j)) is 





(6) 



collision 
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Fig. 3. Real and virtual Feynman diagrams for quark-channel at one-loop order. 



a s N c 
' 2ir 2 

ol s Cf 
2tt 

where 



1-?J (27T 



(x - y) 2 



t/z \ e 



(x-b)l(y-by 



— ik± -r_j_ 



1 



:S^(x x ,y ± ) (7) 



5C 2 )(^,y ± ) = — (U(x ± )UHy ± )), 



N, 



SW(x±,b±,yi 



N, 



2 (Tr[U(x ± )U^b ± )]Tr[U(b ± )U^y ± )}) Y 



In Eq. ([7]) we recognize the collinear divergences associated to the parton distri- 
butions and fragmentation functions, and the rapidity divergence associated the 
dipole gluon distribution. Thus, the collinear divergence is absorbed into the renor- 
malization of the quark distribution [q(x, fi)] 1 ~ loop and fragmentation functions 
[D h / q (z, /x)] 1_loop reproducing the DGLAP evolution equation for each of the dis- 
tribution function. While the rapidity divergence can be absorbed into the renor- 
malization of the dipole gluon distribution [S^(x±, j/±)] 1_loop reproducing in this 
way the Balitsky-BK evolution equation. The factorized formula for the one-loop 
cross section is 

' / '"" ' ~"~ K = f^-^x[ q (x,^-^[D h/q ( Z ,^- 

J Z X 



-loop 



dyd 2 p± 



d 2 x±d 2 yj 



(o) 

2qq 



a »l/(l)' 
2tT 2 «« 



(27T) 



(8) 



Equation © represents the Q CD fa ctorization for hard processes in the satu- 
ration formalism at one loop levef™. The explicit expression for the hard coef- 
ficients can be found in Ref. [TH1 HH1 The factorization formula for inclusive hadron 
production in proton- nucleus collisions is diagrammatically represented in Fig. ([3]). 

It is plausible to think that the QCD factorization formula for inclusive hadron 
production in pA collision obtained at one-loop level holds at any order. In this case 
we can formally write down the expression for the NNLO cross-section as 
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Fig. 4. QCD factorized cross section for inclusive hadron production in pA collisions. On the 
left-hand-side of the picture the black dashed line represents factorization in rapidity at point Y. 
On the right-hand-side of the picture the dashed blue line represents a Wilson line at rapidity Y . 
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We note that at NNLO the parton distributions and fragmentation functions 
follow the NLO DGLAP evolution equation. While the dipole gluon distributions 
follow the Balitsky-JIMWLK evolution equation. In particular we notice that at 
NNLO the distribution follow the NLO B-JIMWLK, S^ 4 ) follow the LO B- 
JIMWLK, while at this order a new operator, (six- Wilson line operators with 
arbitrary white arrangements of color indices), appears and has no evolution. The 
structure of the S^ n ' operators is a hierarchy of evolution equations i.e. the Balitsky- 
hierarchy. 



4. Conclusions 

The main results we have presented here are the NLO fcx-factorization formula 
for DIS and the QCD factorized formula for inclusive hadron production in pA 
collision. We noticed that the kr factorization formula for DIS, Eq. (|6]), is a natural 
consequence of the rapidity factorization in the 2-gluon approximation within the 
OPE at high energy with composite Wilson line operator. Result (O is an analytic 
expression of the NLO impact factor in momentum space. Before, the NLO impact 
factor was available only as a combination of numerical and analytical expressions^! 

In section ([3]) we have presented the QCD factorization result for the inclu- 
sive hadron production in pA collisions. We have shown that the cross section can 
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be written in a factorized form where the parton distributions and fragmentation 
functions follow the DGLAP evolution equation, while the operators S^ n \ due to 
the multiple interactions of the emitted partons with the target, evolve with the 
B-JIMWLK evolution equation. At the end we have also presented a formal expres- 
sion of what would look like the cross section for inclusive hadron production in pA 
collision at two-loop level. 

The author is grateful to the organizers of the QCD evolution 2012 workshop 
for warm hospitality and financial support. 
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